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O ' ABSTRACT: We study a new class of matrix models, formulated on a lattice. On each 

^ , 

0> . site are A^ states with random energies governed by a Gaussian random matrix Hamiltonian. 

d ' The states on different sites are coupled randomly. We calculate the density of and correlation 



between the eigenvalues of the total Hamiltonian in the large N limit. We find that this 
correlation exhibits the same type of universal behavior we discovered recently. Several 



a 

O ■ derivations of this result are given. This class of random matrices allows us to model the 

transition between the "localized" and "extended" regimes within the limited context of 
K^ , random matrix theory. 



We formulate and study a class of matrix models defined in the simplest version by the 
random Hamiltonian 

-C^i I) 

taken from the Gaussian distribution 

p^H) = Le-NtrikimlHf+mlHD+APV^V^^ (2) 

Zj 

Here li\ 2, and V denote A by A matrices, and H a 2 A by 2 A matrix, with H and H\ 2 
hermitean. The parameters m^ 2 ^-nd M are both treated as order A in the large A limit. 



A Hamiltonian of the form in (1) with V equal to zero would naturally arise when a 
conserved symmetry prevents two sectors in the Hilbert space of states to mix. The off- 
diagonal coupling V would then represent symmetry-breaking effects. 

It is natural to generalize our model immediately to contain C sectors which mix due 
to some symmetry-breaking effects. The Hamiltonian /J in (1) is generalized from a 2 by 2 
block matrix to a C by C block matrix. Each block is an A^ by A^ matrix, which we will 
denote by H^^g where a, (3, ... run from 1 to C. The distribution in (2) is generalized to 

Zj 

Here M. represents a C by C real symmetric matrix. The distribution considered in (2) 
corresponds to the case C = 2 with M\y = ml, M22 = "^2' ^^^ -^12 ^ -^21 ^ ^^• 

It is natural also to think of the C sectors as C sites on a lattice. On site a live A^ 
states with energies determined by Haa- The states on sites a and /3 are coupled by H^p. 
The matrix Ai determines the "connectivity" of the lattice. The standard matrix model 
corresponds to the case of a lattice with one site. With this lattice interpretation, this model 
is essentially the same as the model proposed and studied by Wegner some years ago. Our 
results overlap those of Wegner; nevertheless it may be fruitful to approach the model from 
a somewhat different point of view and within the context of our recent study of universal 
correlation in matrix models. '' 

This class of models does not appear to be analyzable by the standard method of orthog- 
onal polynomials. On the other hands, the density and correlation in this class of models 
can be readily obtained using the diagrammatic approach recently developed by us. The 
counting and summing of diagrams provides a somewhat amusing exercise in perturbative 
field theory. 

We note in passing how how we are led to consider such a model. In localization theory 
it is well known that the distribution of spacing between neighboring energy levels changes 
character as one varies the energy from the localized to the extended regime. The physics 
behind this change is clear. Extended wave functions overlap in space, and any perturbation 
would have a significant matrix element between two wave functions close together in energy. 
The familiar phenomenon of level repulsion as described by second order pertubation theory 
tells us that the likelihood of having neighboring energy levels separated by an energy s 
becomes vanishingly small as s goes to zero. This simple argument, essentially due to 
Wigner and to Landau and Smorodinsky, is known as the Wigner surmise. We expect 
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the probability p{s) of a spacing s between neighboring energy levels to rise from zero for 
s = 0, reach a peak, and then to decrease rapidly . The curve traced out by p{s) is sometimes 
referred to informally as "Wigner's head." On the other hand, localized wave functions have 
little overlap with each other, and thus do not repel each other significantly. In the localized 
regime p{s) simply decreases from some finite value as s increases from zero, following a 
Poisson distribution. The curve traced out by p{s) is known informally as "Poisson's tail." 

Recently, there have been some attempts ' to incorporate this transition from the 
"Wigner regime" to the "Poisson regime" in the context of the theory of random matrices. '' 
In matrix models the concept of space does not appear of course; nevertheless one can ask 
whether the distribution of eigenvalues, as manifested in the density of, and the correlation 
between, eigenvalues may exhibit the type of transition described above for localization the- 
ory. In particular, Moshe, Neuberger, and Shapiro proposed a new class of matrix models 
with a parameter they called b such that as b varies the distribution of eigenvalues goes from 
a Wigner to a quasi-Poisson distribution. We found their model somewhat unsatisfactory 
in that the parameter b has to be order A^ larger than the other parameters in the model, 
where A^ as usual in discussions of matrix models denotes the size of the matrices. This 
actually reflects the universality of the Wigner distribution, that is, roughly speaking, its 
tendency to resist change. We feel that the model studied here does not suffer from this 
difficulty. 

Take the C = 2 case as an example. We note that at the special value M = m^ = 
,2 



1712 ~ '^^^eff (which we will refer as the Wigner point), the distribution (2) collapses to 



PiH) = le-2^^^^™^//^^ (4) 

Zj 

the standard Gaussian distribution for IN by IN matrices and we should recover the results 
of Ref. 4 as a check. For M = oo, the sectors 1 and 2 clearly decouples, and the spacing 
distribution is trivially Poissonian in the sense described above. We will refer to this as the 
Poisson point. For M < oo, the off-diagonal perturbation V generates level repulsion, and 
the spacing distribution should be Wignerian. 

We now calculate the density of eigenvalues and the correlation between eigenvalues using 
a diagrammatic method we developed recently. As usual, these quantities can be obtained 
from the Green's functions 

G(z,u.U^(^^tr^^tr^)^ (6) 
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where {0(H)) = J dH 0{H)P{H) and the subscript c indicates the connected Green's func- 
tion. The density of eigenvalues is then given by 

pifi) = (^^tr6ifi -H)'^ = -hmGifi + le) (7) 

and the correlation between eigenvalues, by 

Pc(/x, u) = (^^tr6{fi - H)^tr5{u - H) 

= (-l/47r2)(G'c(++) + Gc{—) - GcX+-) - Gc{-+)) 
with the obvious notation Gc{':^, ±) = Gcl/^ ± «e, z/ ± i5) (signs uncorrelated) . 

The Feynman diagram expansion corresponds to an expansion in powers of 1/ z and 1/w. 
Let us illustrate with the correlation function and expand 

OO OO -. 

{CNfG,{zM =J2Y. -,,^^-^,{trH-trH-), (9) 

The calculation proceeds along the same line as in Ref. 4 suitably generalized to include the 
block indices a, 13,... Diagrammatically, we may borrow the terminology of large N QCD 
and describe the expression for Gc{z,w) as two separate quark loops, of type z and type w 
respectively, interacting by emitting and absorbing gluons. Now each gluon line also carries 
two block indices afi. Some readers may find it natural to think of the C block indices as 
representing a second "color," or "technicolor" index. We take the limit N large for fixed 
G. With the Gaussian distribution in (3) we can readily "Wick-contract" (9). The gluon 
propagator is given by 

{Ha[3ijH^ekl) = J^hl^jk^ae^^p^^afi (10) 

where we have defined the C by C real symmetric matrix a by a^a = 1/M. q. (Note that 
the matrix a is not the inverse of the matrix M. .) Thus, the gluon is represented by a 
double line while a quark is represented by a single line. This convention greatly facilitates 
counting the powers of A^ as explained in Ref. 4. 

Let us begin by calculating the propagator, which is easily seen diagrammatically to be 
diagonal in block indices. Define ga by 

'^aP^ij =< ( ^_ jj )af3,ij >= ^a(3^ij9a (H) 

In the large A^ limit, the sum of the leading planar diagrams ("generalized rainbow" dia- 
grams) is determined by the self-consistent equation 

9a = ^^^ (12) 

^ - Z^/3 ^ap9i3 
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For general Ai this represents a system of coupled non-linear equations to be solved for ga- 
it may be interesting to study the solution numerically for arbitrary A4 . 

We are now ready to calculate Gc{z,w). As in Ref. 4 we begin by ignoring in (9) 
contractions within the same trace (in which case m and n are required to be equal). In 
the large N limit, the dominant graphs are given essentially by "ladder graphs" (with one 
crossing). We now have a combinatorial problem which is most easily solved by distorting 
the ladder graph in Fig. (la) to the "wheel" graph in Fig. (lb). (Here we admit distortions 
as long as they preserve the combinatorial factor.) We see that each "spoke" of the wheel is 
associated with a factor of o"Q,fl. Thus, the graphs with n loops should be multiplied by the 
factor tro"". We obtain 

iCNfGciz,w) = f^^^^^ (13) 

We next include Wick-contractions within the same trace in {trH^trH'"'). Graphically 
these contractions correspond to decorating the ladder graphs by vertex and self energy 
corrections. First, we have to correct the upper part of the graph in Fig. (la) with vertex 
corrections, as shown in Fig. (2), and similarly for the lower part of the graph. These vertex 
corrections lead to a rather complicated expression. For example, the factor tra^ in (13) for 
n = 3 is to be replaced by the sum of terms, one of which has the form 

^-^ ^ 1 — a z'^ ^ ^-^ 1 — a w^ 

apuj (3 1 

and another with the factor ^ ( -._ , 2 )07 i^i the expression above replaced by ^^ ( ,_ , 2 )cj7- 
This clearly leads to a rather unwieldy expression into which we still have to put in the self 
energy corrections. 

Clearly we should not expect to obtain a simple and closed expression for arbitrary a: 
every sector or every site would be unique in its local properties. Instead, we made the 
rather mild and physically reasonable assumption of homogeneity by supposing that 

y ^(JqP = independent of a. (15) 

Every site (or sector, if the reader prefers) is treated on the same footing. In particular, this 
assumption holds if the lattice is translation invariant. 

Quite remarkably, we now notice that this model may in fact be solved on any lattice in 
any dimension as long as the condition (15) holds. Consider the eigenvalue decomposition 
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the real symmetric matrix cr, which we may identify as a "hopping matrix" , 

a = Y,\k>ek<k\ (16) 

k 

where k takes on C values. Denote the eigenvalues by e^. The condition (15) merely states 
that of the eigenvectors of a there is one (which we labeled as |A; = >) with components 
all equal to Xj \fC . 



We find immediately that in (12) g^ is independent of a and equal to 



G{z) = ^ -; ^^Q (17) 

2eo 

This has the same form as the propagator in the simple one matrix model and gives imme- 
diately Wigner's semi-circle law for the density of eigenvalues 
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p(/^) = ^V^eo-/i2 (18) 

We now proceed to the two-point Green's function. The vertex correction factor we 
encountered before now simplifies drastically 

E( i \ 2 W = ( . \ 2 ) (19) 

to a factor independent of a. Putting in vertex corrections just amounts to multiplying by 
this factor and a corresponding one with z replaced by w. The trace and sum in (13) can be 
performed immediately to give 

{CNfGeiz, w) = J2 , '' ,2 - (20) 



After putting in vertex and self-energy corrections, we find 

G'^iz) \ ( G^[w) 

k 



VI - eQG^{z)J\l- eQG^{w)J V^ (1 - ekG{z)G{w)yJ 



Taking the appropriate absorptive parts we can calculate the connected correlation func- 
tion pcifJ', J^)- As in our earlier work we find it natural to introduce angular variables defined 
by 

sin^ = -^ (22) 

2eo 
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and 

V 

smv9 = -— . (23) 

2eo 

The angles Q and (/? vary from — 7r/2 to +7r/2 over the range of the eigenvalues. The connected 

correlation function is given by 

1 



Pc{^. v) 



IGvr^A^^C^eo cos Q cos Lp 
^ ^ rfc + cos(g + y) rfc-cos(g-v?) -| (^4) 

^ ^1 a + rfcCOs(^ + 99))2 (l_rfcCos(^-v5))2, 



where we have defined 

ru = |2^^ (25) 

^0 + ^fc 

The expressions (21) and (24) represent the central results of this paper. It is worth empha- 
sizing that our results hold for any lattice with any connectivity as long as the homogeneity 
condition (15) is satisfied. 

As a simple special case, consider the simple case in which the matrix A^ is such that all 
its diagonal elements are equal to m and all its off-diagonal elements are equal to M . Every 
site on the lattice couples to every other site. The matrix a has one eigenvalue eg = l/'o^f 
and C — 1 degenerate eigenvalues equal to l/""^^- Thus we have 



where 



D ia -) - ml f 1 1 




'"-^'^''^' 16n^N^C^cos9cos<^il + cos{e + ip) 1 - cos(^ - 


-v) 


r + cos(^ + V9) ^_cos(^-v^ ) 




^ ^■(l + rcos(^ + V?))2 ' ^"^ ^■(l-rcos(^-V^))2J 




2m1ml 
mj + mf 





(26) 



(27) 



Note that |r| < 1 and thus only the second term in pc{^.,v) is singular when fi = u. As 
we vary M , the correlation moves from a Poissonian ("localized") regime to a Wignerian 
("extended") regime. In particular, at the Wigner point, r = and we recover trivially the 
result for the one matrix model first obtained in Ref. 2. At the Poisson point, r = 1 we 
have a short distance (i.e., when 6 and (/? approach each other) singularity with a residue 
(C — 1 + 1)/C . (Recall from Ref. 2 that this short distance singularity occurs because we are 
dealing with the "smoothed" correlation. The exact correlation does not have this singularity. 
Nevertheless, this singularity tells us about the short distance correlation between the density 
of eigenvalues.) Thus, the short distance singularity is reduced from the standard C = 1 
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case and is completely suppressed in the large C limit, namely in the limit of many sectors 
forbidden to mix by many conservation laws. More generally, for any r, it is sufficient to 
take C large for this short distance singularity to be suppressed. 

For an arbitrary lattice, we simply find the eigenvalues of the hopping matrix a. For 
example, for a hypercubic lattice we have 

1 2 . 

efc = ^ + ]^E^°^^« (28) 

a 

where the sum over a runs overs the dimension D of the lattice. In the expressions (21) and 
(24) the sum over k runs as usual over the Brillouin zone ka = 27rja/La,ja = 0, 1, 2, ...La — 1 
with LiL2...L]j = C. Note that to have nearest neighbor coupling the appropriate entries 
in the matrix M. have to be taken to infinity. This simply means that some of the matrices 
Hfj^fj in (3) are to be set equal to zero. It may be interesting to study our result for a variety 
of lattices. Also, one may consider the limit C ^ oo in which case the sum over k is replaced 
by an integral. While our results are obtained in the context of random matrix theory, they 
may conceivably be of relevance to some condensed matter systems. 

This correlation function, while it is not quite the same as the correlation functions 
we ' ' and others ' have found previously in a number of different situations, has the 
same general structure as these other correlation functions. This in itself is a statement of 
universality. Next, we may make a remark about going beyond Gaussian distributions. In 
our previous work ' we identified two classes of matrix models: the Wigner class and the 
trace class. Briefiy, in the Wigner class the individual matrix element obeys a probability 
distribution. (In the present context, the distributions for different H^^g may in general be 
different.) As was explained in Ref. 3, universality, that is, the independence of the density 
and the correlation on these probability distributions, can be proved almost immediately 
using our diagrammatic method. In contrast, in the trace class, each of the H^^g obeys some 
probability distribution of the form 

P{Hap) <x e-'V"-l~''(«»») (29) 

(no sum over repeated indices) where V^^ denotes a polynomial, different for different ap. 
In this case, we do not know how to prove universality, that is, the independence of the 
correlation function pc (when expressed in terms of scaling variables) on V*^^ . It would be 
interesting to see numerically whether universality indeed holds. 

We would now like to present two other derivations of our central result. Let us for 
notational simplicity consider the case studied in Ref. 4, that of one Gaussian-random matrix. 



which we will denote by H . In other words, we first consider H for C = 1, a lattice with 
one point. We will first recover the results of Ref. 4 and then generalize these results to an 
arbitrary lattice. Consider the "scattering amplitude" 



A^(((^)K^)™) - 51G(.)CGH) 



^ 6l,6rA + 616^ B 

The two "scalar" functions A and B depend on z and w of course. The Green's function is 
defined as usual by 



we now calculate (repeated indices summed) 

ril = N'^A + NB = N^[^^^^^^^^-G{z)G{w)] (32) 



This equation tells us that A is of order A^ while B is of order 1/A^. Furthermore, the 
connected two-point Green's function is given by 

N'^Gciz, w) = (tr^-—tr^-—) =A + NB (33) 

\ z-H w-H/^ 

Thus, if we know G and Gc we know A and B and hence the full scattering amplitude r|"^. 
For our lattice problem we have to generalize the definition above to 






(34) 



All is as before except that now the scattering functions A and B can depend on the site 
labels a, /9, 7 and e or technicolor indices because the probability distribution is not unitary 
invariant in these indices. Note that nevertheless the kronecker delta structure in (34) holds. 
This can be seen graphically since a line carrying color and technicolor going in must come 
out carrying the same color and technicolor. 

Now we note the remarkable fact that if we are interested in questions of localization, 
namely what happens when the site a is equal to the site e, and far away from the site /5=7, 
we don't need to know anything about B. But A is determined by the Green's function G{z) 
completely! 
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Recalling that for the Gaussian and single matrix case (see (12) and (17)) that G{z) is 
the solution of the quadratic equation 



G + i (35) 



we obtain immediately from (32) that 



^'-•"'= l-G(.)GW <=*«' 

One can see readily that for the case of a lattice of matrices considered here we only have to 
modify this result to 

Note that as anticipated A depends on a and j3. 

This result appears to agree with Wegner's result . We would like to emphasize the 
simplicity and ease with which this result can be obtained. For the sake of completeness, let 
us recall that the standard localization formalism instructs us to set 2; = i? + iu and w = 
E—iuj and then let the distance r = |a— /3| tend to infinity. Since G{E +iuj)G{E —ioo) ^ 1 as 
cij ^ and e^^ — > 1 — 0{k) as A; ^ for a simple rectilinear lattice, we see that A -^ 1/r , 
as is well known from the work of Wegner. 

As a check, we may also easily calculate A directly. From the structure of the indices we 
see that A is given by summing the ladder graphs in Fig (3). Note that vertex corrections 
are suppressed. We merely have to dress the quark propagators. We thus have 

A(z,^) = — 5^^(0^ (38) 

n—1 

After summing and dressing the quark propagators by replacing 1/ z by G{z) and 1/w by 
G{w) we obtain (37) immediately. 

In this approach, in order to obtain our correlation function, we must calculate B. Thus, 
there is conservation of labor after all. While it is somewhat more involved to calculate 
B than to calculate A, it is still easy enough. The relevant graphs are given in Fig. (4). 
Curiously, because of the structure of indices only graphs with vertex corrections contribute 
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to B. In particular, the simple one-gluon exchange graph does not contribute. We obtain 
(again for notational simplicity we do the C = 1 case here) 



OO OO _. OO -. 



^^(-'-) = (EEt2=-i)E 



z'^Pw'^1 ' ^ (zw)"+2 

^ {1-G^{z)){l-G^{w)) 1 - G{z)G{w) 

after dressing the quark propagators. Note the (—1) inside the parenthesis. We leave it to 
the reader to show that when this somewhat strange looking expression is combined with 
A according to (33) we recover our result ' for the correlation function Gc{z,w). It is a 
straightforward exercise to put in the a matrix to recover (21). 

We conclude by giving yet another derivation of our central result. Note that (6) may 
be written as 

Again, for notational simplicity, let us consider the G = 1 case first. Expanding the loga- 
rithms, we find 

9 9^^ 

n=l m=l 

This is represented by the RwheelS graph of Fig (lb). We see immediately that vertex 
corrections are suppressed in the large N limit. Thus, we may set n = m in (41) and since 
(trH^trH^) = n we can immediately evaluate the sum to obtain 

N^G{z, w)c = -^-^ log(l - — ) (42) 

oz ow zw 

Dressing the quark propagators, we obtain immediately 

N^Giz, w)c = -§;^ log(l - Giz)Giw)) (43) 

This derivation is simpler than that given in Ref. 4. We also obtain a more compact form 
as given in Ref. 18. It is now simple to go to the case of an arbitrary lattice. We now have 
{trH'^'trH'"') = ntra^ and thus we find 

N^G{z, w)c = -^^trlog(l - aG{z)G{w)) = -J^ ^^^^^(l " ^kG{z)G{w)) (44) 

k 
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Gc(z,w) = ^TT-V V -^^( —trH^'^trH'^) (41) 



We leave it to the reader to check that we have indeed recovered our previous result (21). 

Finally, we are tempted to indulge in a speculation. It has long been known that matrix 
models are intimately related to the Calogero-Sutherland model. Specifically, the Calogero- 
Sutherland model for a certain coupling constant A equal to 1/2, 1, and 2 correspond re- 
spectively to matrix models with real symmetric, hermitean, and quarternionic matrices. 
Recently, Ha was able to solve the Calogero-Sutherland model for A equal to an arbitary 
rational number. The question naturally arises whether or not for these values of the cou- 
pling constant the Calogero-Sutherland model corresponds to matrix models. Is it possible 
that the class of matrix models considered here may be the correspondents? And if not, it 
would be interesting to ask whether the matrix models considered here corresponds to some 
generalized Calogero-Sutherland models. 
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FIGURE CAPTIONS 

Fig (la) A crossed ladder graph that contributes to leading order in the large A^ limit; 
(lb) As far as the combinatorics is concerned, graphs of the type drawn in (la) are equivalent 
to the graphs of the type drawn in (lb). 

Fig (2) Vertex corrections to the upper part of the graph in (la) 

Fig (3) Graphs contributing to A in the large A^ limit. 

Fig (4) Some graphs contributing to B; note that only graphs with vertex corrections 
are to be included. 
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